STEIN'S METHOD, HEAT KERNEL, AND TRACES OF 
POWERS OF ELEMENTS OF COMPACT LIE GROUPS 

JASON FULMAN 

Abstract. Combining Stein's method with heat kernel techniques, we 
show that the trace of the jth power of an element of U(n, C), USp(n, C) 
or SO(n, R) has a normal limit with error term of order j/n. In con- 
trast to previous works, here j may be growing with n. The technique 
^-^, should prove useful in the study of the value distribution of approximate 

eigenfunctions of Laplacians. 
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1. Introduction 



There is a large literature on the traces of powers of random elements 
of compact Lie groups. One of the earliest results is due to Diaconis and 
Shahshahani [5J. Using the method of moments, they show that if M is ran- 
dom from the Haar measure of the unitary group U(n, C), and Z = X + iY is 
a standard complex normal with X and Y independent, mean and variance 
2 normal variables, then for j = 1, 2, • • • , the traces Tr(M^) are indepen- 
dent and distributed as \fjZ asymptotically as n — > oo. They give similar 
results for the orthogonal group 0(n,~R) and the group of unitary symplectic 
q | matrices USp(2n,C). The moment computations of [5j use representation 

C*~) ■ theory. It is worth noting that there are other approaches to their moment 



computations: [19j] uses a version of integration by parts, [ll[ uses the com- 
binatorics of cumulant expansions, and [3( uses an extended Wick calculus. 
We mention that traces of powers of random matrices have been studied for 
other matrix ensembles too ( [2(| . |7|| . [28l] V 

Concerning the error in the normal approximation, Diaconis conjectured 
that for fixed j, it decreases exponentially or even subexponentially in n. 
/\ ' In an ingenious paper (which is quite technical and seems tricky to apply 

to other settings), Stein [30[] uses an iterative version of "Stein's method" 
to show that for j fixed, Tr(M J ) on Q(n,M) is asymptotically normal with 
error 0{n~ r ) for any fixed r. Johansson [12] proved Diaconis' conjecture for 
classical compact Lie groups using Toeplitz determinants and a very detailed 
analysis of characteristic functions. Duits and Johansson [6( allow j to grow 
with n in the unitary case, but do not obtain error terms. We also note that 
in the unitary case when j > n, the situation is not so interesting, since by 
work of Rains [2j|, the eigenvalues of M J are simply n independent points 
from the unit circle (and he proves analogous results for other compact Lie 
groups) . 



Date: May 4, 2010. 

Key words and phrases. Random matrix, Stein's method, heat kernel. 

Dedicated to Thuy Le, on the occasion of our tenth anniversary. 

1 



2 JASON FULMAN 

The current paper studies the distribution of Tr(M 3 ) using Stein's method 
and heat kernel techniques. This is a follow-up work to the paper [9[], 
which used Stein's method and character theory to study the distribution of 
x(M), where x 1S the character of an irreducible representation; the func- 
tions Tr(M J ) are not irreducible characters for j > 1, so do not fit into the 
framework of [9(. It should also be mentioned that the heat kernel is a truly 
remarkable tool appearing in many parts of mathematics (see the article 
[13j ] for a spirited defense of this statement with many references), and we 
suspect that the blending of heat kernel techniques with Stein's method will 
be useful for other problems. 

In earlier work, Meckes [17j], used Stein's method to study eigenfunctions 



of the Laplacian (a topic of interest in quantum chaos and arithmetic 27], 
among other places). We note two differences with her work. First, she 
uses geodesic flows and Liouville measure instead of heat kernels. Second, 
her infinitesimal version of Stein's method [171 ]. [181 ] uses an exchangeable 
pair of random variables (W,W e ) with the conditional expectation E[W e — 
VF|VF] divided by e 2 approximately proportional to W as e — > 0. In the 
current paper the natural condition is that E[W e — VK|VK] divided by e is 
approximately proportional to W as e — > 0. 

We do use some moment computations from [5j, but as is typical with 
Stein's method, only a few low order moments are needed. It should also 
be mentioned that the constants in our error terms can be made completely 
explicit (for instance in the unitary case we prove a bound of — -^), but we 
do not work out the other constants as the bookkeeping is tedious and the 
true convergence rate is likely to be of a sharper order. As to future work, 
we note that more general linear combinations of traces of powers do satisfy 
central limit theorems (see [J], [121 ]. [291 ] for precise conditions); obtaining 
good error terms by our techniques (or other methods) may be quite tricky 
and is an important problem. 

The organization of this paper is as follows. ISection "21 gives background 
on both Stein's method and the heat kernel. ISection "31 treats the orthogo- 
nal groups. fSection 4l treats the symplectic groups, and ISection 51 treats the 
unitary groups. 

2. Stein's method and the heat kernel 
In this section we briefly review Stein's method for normal approximation, 



using the method of exchangeable pairs 3l| . One can also use couplings to 
prove normal approximations by Stein's method (see [22] for a survey), but 
the exchangeable pairs approach is effective for our purposes. For a survey 
discussing both exchangeable pairs and couplings, the paper [23J can be 
consulted. 

Two random variables W, W on a state space X are called exchangeable 
if the distribution of (W, W) is the same as the distribution of (W , W). As 
is typical in probability theory, let E(A|_B) denote the expected value of A 



given B. The following result of Rinott and Rotar 24j| uses an exchangeable 



pair (W, W) to prove a central limit theorem for W . 
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Theorem 2.1. ([21]) Let (W,W) be an exchangeable 'pair of real random 
variables such thatE(W) = 0,E(W 2 ) = 1 andE(W'\W) = (l-a)W + R{W) 
with < a < 1. Then for all real xq, 



P(W < 
p. 


1 rxo 

\/27r 7_oo 


x 2 

' 2 dx 


-^Var 
a 


(E[(W> W) 2 \W\) + 19 



£ -JVar(E[(W - W) 2 \W\) + 19 V E ^ I + 6 A /iE|H" - Wl 3 . 
a a V a 

In practice, it can be quite challenging to construct exchangeable pairs 
satisfying the hypotheses of Theorem 12.11 an d such that the error terms are 
tractable and small. 

Lemma 12.21 is a known inequality (already used in the monograph [31]) 
and useful because often the right hand sides are easier to compute or bound 
than the left hand sides. We include the short proof. Here M is an element 
of the state space X (in this paper X is a compact Lie group and M a matrix 
inX). 

Lemma 2.2. (1) Var(E[(W - W) 2 \W}) < Var(E[(W - W) 2 \M}). 

(2) With notation as in Theorem\2j\ letting E(W'\M) = (1 - a)W + 
R(M), one has that E(R{W) 2 ) < E(R(M) 2 ). 

Proof. Jensen's inequality states that if g is a convex function, and Z a ran- 
dom variable, then g(E(Z)) < E{g(Z)). There is also a conditional version of 
Jensen's inequality (Section 4.1 of [8]) which states that for any a subalgebra 
F of the a- algebra of all subsets of X, 

E(g(E(Z\F))) < E(g(Z)). 

Part 1 now follows by setting g(t) =t 2 , Z = E((W' - W) 2 \M), and letting 
F be the a-algebra generated by the level sets of W . Part 2 follows by setting 
g(t) = t 2 , Z = R(M), and letting F be the cr-algebra generated by the level 
sets (AW. □ 

To construct an exchangeable pair to be used in our applications, we use 
the heat kernel on G. See [iCJ], [25(, 26|], [16|]) for a detailed discussion of 



heat kernels on compact Lie groups, including all of the pro per ties stated 
in the remainder of this section. The papers [LJ], [l|], [20|, [l5| illustrate 
combinatorial uses of heat kernels on compact Lie groups. 

The heat kernel on G is denned by setting for x,y £ G and t > 0, 

(1) K(t, x,y) = Y, e- Xnt Mx)<Pn(y), 

n>0 

where the X n are the eigenvalues of the Laplacian repeated according to 
multiplicity, and the (j> n are an orthonormal basis of eigenfunctions of L 2 (G); 
these can be taken to be the irreducible characters of G. 

We use the following properties of the heat kernel. Here A denotes the 
Laplacian of G, and e is defined as I + £A + t 2 -y- + • • • . 

Lemma 2.3. Let G be a compact Lie group, x,y € G, and t > 0. 

(1) K(t,x,y) converges and is non-negative for all x,y,t. 
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(2) J G K(t, x, y)dy = 1, where the integration is with respect to Haar 
measure of G. 

(3) e tA (p(x) = f yeG K(t,x,y)4>(y)dy for smooth <j). 

The symmetry in x and y of K (t, x, y) shows that the heat kernel is a 
reversible Markov process with respect to the Haar measure of G. Thus, 
given a function / on G, one can construct an exchangeable pair (W, W) by 
letting W = f{M) where M is chosen from Haar measure, and W = f(M'), 
where M' is obtained by moving time t from M via the heat kernel. 

3. The orthogonal group 

If A is an integer partition (possibly with negative parts) and mj denotes 
the multiplicity of part j in A, we define p\{M) = JT ■ Tr(M 3 ) m i . For exam- 
ple, P5,3,3(M) = Tr(M 5 )Tr(M 3 ) 2 . Typically we suppress the M and use the 
notation p\. We let W = ^4* if j is odd and let W = "^-p- if j is even. Note 
that since the eigenvalues of M are roots of unity and come in conjugate 
pairs, pj = p-j is real. The main result of this section is a central limit 
theorem for W with error term 0(j/n). 

The following moment computation of [11] (analogous to that of [5] for 
the full orthogonal group) will be helpful. In fact as the reader will see, in 
the applications of Lemma 13.11 we only use fourth moments and lower. 

Lemma 3.1. Let M be Haar distributed on SO(n,M). Let (ai,02,--- , Ofe) 
be a vector of non-negative integers. Let Z±, ■ • ■ ,Zk be independent standard 
normal random variables. Let r]j be 1 if j is even and otherwise. Then if 

n-l> Yll=i <H, 



E 



Y[Tr(M j ) a i 



k k 



Here 



. „ . . , , , N f if a is odd 

if j is odd, gjla) = <.„/•>, n w oX n .„ 
J ' y3K ' \ fl\a- l)(a-3)---l if a is even 

if j is even, 9j {a) = 1 + J^ ( 2 lV( 2fc " X )( 2fc - 3) ■ ■ ■ 1- 

fe>i ^ ' 

Rains [20(] (see also [lj]) determined how the Laplacian acts on power sum 
symmetric functions. We need his formula only in the following two cases. 

Lemma 3.2. (1) 

(n — l)j j TT-^ j TT-^ 

&SO(n)Pj = 2 p i ~ 2 ^ Pl ' j ~ l + 2 ^ P2l ' j ' 



&so(n)Pj,j = -{n- l)j PjJ - j 2 p 2j - JPj Y, PiJ-i + 3Pj Yl P2l -i + j 



2 n. 
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We fix t > 0, and motivated bv lSection 2\ define 

f k 
W' = e tA (W) = W + Y^ ry Afc (W). 



k>i 



Lemma 13.31 computes the conditional expectation E[W'|M]. 
Lemma 3.3. 



E[W'\M] = ( 1 - ^ — ^ ) II ' + R(M), 



with 

R{M) = t 

and 
R(M) = t 



Vl 



Y Pu-i + ^r £ 



P21-J 



1<1<3 



1<1<3 



+ 0(t 2 ) j odd, 






+0(t 2 ) j even. 



Proof. Applying part 3 of Lemma 12.31 and part 1 of Lemma [37 

E[W'\M] 
= e tA (W) 



W + t 



2 Pj " T ^ Pl ' j - L + ~2 ^ P2 '" 



i<z<i 



i<z<j 



+ 0(* 2 ), 



and the result follows. 



□ 



Lemma l3~4l computes E[(W' — W) 2 |M], a quantity needed to apply The- 
orem [2J] Many cancelations occur, and a simple formula emerges. 

Lemma 3.4. 

E[(W - W) 2 \M] = tj(n - p 2j ) + 0{t 2 ). 
Proof. Clearly 

E[(W' - W) 2 \M] = E[(W') 2 \M] - 2WE[W'\M] + W 2 . 
Suppose now that j is odd. By part 3 of Lemma 12.31 and part 2 of Lemma 

E[{W') 2 \M] 
= W 2 + -A Pjd + 0(t 2 ) 



W 2 + t 
+0{t 2 ). 



-(n - l)pjj - j>2i - Pj Yl PU-i + Pi Yl P2l ~i + J ' 



n 
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By Lemma E31 -2WE\W'\M\ is equal to 



-2W 2 +t 



(n-l)jW 2 + Pj Y^ Pl,j- 
l<l<3 



Pi Yl P2l -i 

l<Kj 



+ 0{t 2 ). 



Thus 



E[(W') 2 \M\ - 2WE[W'\M] + W 2 = tj{n - p 2j ) + 0{t 2 ), 

as claimed. A very similar calculation shows that the same conclusion holds 
for j even. □ 

Lemma 3.5. Suppose that 4j < n — 1. Then 

Var(E[(W' - W) 2 \M]) = 2ft 2 + 0(t 3 ). 

Proof. By Lemma 13.41 

Var(E[(W' - W) 2 \M\) = j 2 t 2 Var(p 2j ) + 0(t 3 ). 

The result now follows from Lemma 13. 11 □ 

Lemma 3.6. Suppose that 4j < n — 1. Then 

(1) E(W - W) 2 = tj(n - 1) + 0(t 2 ). 

(2) E(V' -W) i = 0{t 2 ). 

Proof. Lemma E31 implies that E(W - W) 2 = E[tj(n-p 2j )] + 0(t 2 ). From 
Lemma 13-H E(p 2 j) = 1, which proves part 1 (only using that 2j < n — 1). 
For part 2, first note that since 

E[(W'-Wf] = E(W 4 )-AE{W 3 W')+QE[W 2 {W') 2 ]-AE[W{W') 3 ]+E[(W')% 

exchangeability of (W, W) gives that 

E{W'-W) A = 2E(W A )-8E(W 3 W') + QE[W 2 (W') 2 ] 

= 2E(W A ) - 8E[W 3 E[W'\M\] + 6E[Ly 2 E[(W') 2 |M]]. 
Supposing that j is odd and using Lemma 13.21 this simplifies to 

2E(W 4 ) - 8E[W 4 ] + 6E[W 4 ] 



+tE 



^n-l^W^ + AW 3 ^] J2 Pid-i-^W 3 yfj Y, P21- 

l<Kj l<l<j 



-tE 



-6(n - l)jW 4 - 6W 2 Pj Y, PU-i + m2 Pi 2 P»-J 

i<z<j i<2<j 

+tE [-6jW 2 p 2j + 6VF 2 jra] + 0(t 2 ). 
By Lemma 13. 11 this simplifies to 

t [I2j(n - 1) - 18j(n - 1) - 6j + 6jn] + 0(t 2 ) = 0(i 2 ), 

as claimed. A very similar calculation gives the same conclusion for j even. 

□ 



Next we bound a quantity appearing in the second term of Theorem 12.11 
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Lemma 3.7. Suppose that 4j < n — 1. Let 



R(M) = t 
and 
R(M) = t 






+ 0(t 2 ) j odd, 



n IT Z> Vl '3~l + IT 1^ PU-3 



l<Kj 



l<Kj 



+0(t ) j even. 



ThenE[R 2 } = 0(t 2 f). 



Proof. Suppose that j is odd. Applying Lemma 13, II and keeping only terms 
with non-0 expectations, one has that 



E[R 2 



t 2 J 



E 



4 E (pu-i) 2 + 4 E (^) 2 - 8 E wj-i 



I odd 



I odd 



I odd 



+ 0(t 3 ) 



t 2 j 



4 e iu - 1 + 1) - 4 E l 



I odd 



I odd 



+ 0(t 3 ) 



= 0(t 2 f). 
The case of j even is proved in a similar way, as can be seen by writing 



R = t 



"2 2" ^ W,i ~' ~2~ ^-> P2/_i 



i<*<? 



i<kj 



+ 0(t 2 ). 



D 



Combining the above calculations leads to the main result of this section. 

Theorem 3.8. Let M be chosen from the Haar measure of SO(n,M). Let 
W(M) = llipi ifj i s dd and W(M) = TliMpzl if j is even . Then 



\W < x ) 



1 



e 2 dx 

2-7T J -oo 



0(j/n). 



Proof. The result is trivial if 4j > n — 1, so assume that 4j < n — 1. We 



Kn-l)j 



and 



apply Theorem 12.11 to the exchangeable pair (W, W') with a 
will take the limit t — > in each term (keeping j,n fixed). By part 1 of 
Lemma [2.2l and Lemma [3.51 the first term is 0(\/~j/n). By part 2 of Lemma 
12.21 and Lemma 13.71 the second term is 0(j/n). By the Cauchy-Schwarz 
inequality and Lemma |3.6[ 



E\W' - W\ 3 < ^E{W - W) 2 E(W - wy = 0{t 3/2 ). 

Thus the third term in Theorem 12.11 tends to as t — > 0, and the result is 
proved. □ 
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with all blocks 



4. The symplectic group 

Let J be the 2n x 2n matrix of the form I T 

n x n. USp(2n, C) is defined as the set of 2n x 2n unitary matrices M with 
complex entries such that MJM 1 = J; it consists of the matrices preserving 
an alternating form. As in ISection "3l we use the notation that p\{M) = 
n • Tr(M J ) m J , and we typically suppress the M and use the notation p\. 

We let W = 2j* if j is odd and let W = Pj ~ if j is even. Since the 
eigenvalues of M are roots of unity and come in conjugate pairs, pj = p-j 
is real valued. The main result of this section is a central limit theorem for 
W with error term 0(j/n). 

The following moment computation is the symplectic analog of Lemma 
13.11 It was proved by [5( under the slightly weaker assumption that n > 
Yli=i a k- As stated, Lemma I4~T1 appears in [ll|, with a later proof in [191 ] . 

Lemma 4.1. Let M be Haar distributed on USp(2n, C). Let (ai, 02, • • • , ak) 
be a vector of non-negative integers. Let Z±,--- ,Zk be independent standard 
normal random variables. Let r/j be 1 if j is even and otherwise. Then if 

2n + l > J2i=i a i> 

k 

E 



Y[Tr(M j ) a > 
3=1 



n(-i)°'- i)a ^,(a J )=n iE (v^^ 



i=i 



Vj, 



Rains [20J (see also [14]]) determined how the Laplacian acts on power sum 
symmetric functions. We need his formula only in the following two cases. 



where the polynomials gj are as in Lemma \3. 1 



(2n + l)j 



- y 

2 *-" 



P21-J 



l<Kj 



Yl pij-i- 

l<Kj 



Lemma 4.2. (1) 

^USp(2n)Pj = ~- 2~~^ Pj 

(2) 

A USp(2n)Pj,j = -(^ + l)jPjj-j 2 p2j-jPj J2 VlX-j-jVj 5Z Pld-l + 2 f 

As in the orthogonal case, we fix t > 0, and define 

t k 



n. 



W' = e tA (W) = W + Y, ry Afc W- 



fc>i 



Lemma 4.3. 



E[W'\M] = ( 1 - t(2w ^ 1)j ) W + R(M), 



with 



R{M) = t 






l<Kj 



l<Kj 



+ 0(t 2 ) j odd, 
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and 
R(M) = t 



(2n + l)Vj VJ ^ VJ ^ 
o T 2^ P2l -i ~ T L. PU- 



+0(t ) j even. 



Proof. Applying part 3 of Lemma 12.31 and part 1 of Lemma 14,2 

E[W'\W] 
= e tA (W) 



W + t 



(2n + l)Vj VI 



-Pj 



Vl 



VJ v^ VJ v^ 



mj-i 



1<KJ 



l<l<j 



+ 0(t 2 ), 



and the result follows. 



□ 



Lemma |4~41 computes E[(W' — PF) 2 |M], a quantity needed to apply Theo- 
rem [2H1 As in the orthogonal case, there are many cancelations, leading to 
a simple formula. 

Lemma 4.4. 

E[(W' - W) 2 \M] = tj (2n- P2j ) + 0(t 2 ). 

Proof. Clearly 

E[(W' - Wf\M] = E[(W') 2 \M] - 2WE[W'\M] + W 2 . 

Suppose that j is odd. By part 3 of Lemma 12.31 and part 2 of Lemma |4.2[ 

E[(W') 2 \M] 

= W 2 + -A PjJ + 0(t 2 ) 



H ' I- / -(2n + l)pjj - jp 2j - Pj ^2 Pv-3 ~ Pi ^2 Wj- 1 + 2 J 

l<l<j l<l<j 

+0(t 2 ). 

By Lemma O -2WTE[W|M] is equal to 



n 



-2W 2 + t 



(2n + l)pjj+Pj Yl Pu-i + Pj £ PiJ-i 

i<i<j i<i<j 



+ 0(t 2 



Thus 



E[(W') 2 \M] - 2WE[W'\M] + W 2 = tj [2n - p 2j ] + 0(t 2 ) 
as needed. A similar computation proves the lemma for j even. 

Lemma 4.5. Suppose that 4j < 2n + 1. Then 

Var(E[(W' - W) 2 \M]) = 2ft 2 + 0{t 3 ). 
Proof. By Lemma 14.41 

Var(E[(W' - W) 2 \M\) = j 2 t 2 Var(p 2j ) + 0{t 3 ). 
The result now follows from Lemma [4.11 



□ 



□ 
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Lemma 4.6. Suppose that Aj < 2n + 1. 

(1) E(W - W) 2 = tj(2n + 1) + 0{t 2 ). 

(2) E{W -Wf = 0{t 2 ). 

Proof. Lemma 031 implies that E(W - VK) 2 = E [tf (2n - p 2 j)} + 0(t 2 ). 
From Lemma 14.11 E(p 2 j) = — 1, which proves part 1 (even assuming that 
2j<2n + l). 

For part 2, first note that since 

E[(W'-W) A ] = E(W A )-AE{W 3 W')+6E[W 2 {W') 2 ]-4E[W{W') 3 ]+E[{W')% 

exchangeability of (W, W) gives that 

E(W' -W) A = 2E(W A )-8E(W 3 W) + 6E[W 2 (W') 2 ] 

= 2E(W A ) - 8E[VF 3 E[M/'|M]] + 6E[W 2 E[(W') 2 \M}]. 

Suppose j is odd. Using Lemma r4.3l and part 2 of Lemma l4.2l this becomes 

2E{W i ) - 8E(W 4 ) + 6E(VF 4 ) 

+£E 



+tE 



4(2n + l)jW 4 + AW 3 ^j J2 PU-i + ^W 3 ^j ^ p 2l ^ 

l<Kj l<Kj 



-6(2n + l)jW 4 - 6jW 2 p 2j - 6W 2 Pj ^ p y _, 

1<1<3 



+t 



-6W 2 Pj Y, P2i-j + ^W 2 j 



i 'i 



+ 0(t 2 



By Lemma l4.lt this simplifies to 

t [12j(2n + 1) - 18j(2n + 1) + 6j + 12jn] + 0(t 2 ) = 0(t 2 ), 
as claimed. A similar calculation gives the same result for j even. 
Lemma 4.7. Suppose that 4j < 2n + 1. Let 



□ 



R(M) = t 



VJ 



VI 



VJ v^ VJ v^ 



J'-» 



i<z<j 



1<1<3 



+ 0(t 2 ) j odd, 



and 



R(M) = t 



(2n + l)^ V7 



E 

l<Kj 



P21-J 



VI 



E 



Pl,j-l 



+0{t 2 ) j even. 



Then E[R 2 ] = 0{t 2 j A ) 
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Proof. Suppose that j is odd. Applying Lemma 14, II and keeping only terms 
with non-0 contribution, one has that 



E[R 2 



t 2 j 



E 



4 ^(w,^) 2 + 4^(w) 2 + 8 X> M ,;-* 



I odd 



I odd 



I odd 



o(t 3 ) 



t 2 j 



4 Y, Z(j - Z + 1) - 4 Y I 



I odd 



I odd 



+ 0(t 3 ) 



2„-4> 



o(?j 



The case of j even is proved by a similar argument, after writing 

+ 0(t 2 ). 



R = t 



V7 VJ sr 

"2 T L P«-i 



i<kj 

W/2 



2 



i<z<j 



P/J-Z 



D 



Theorem 4.8. Let M be chosen from the Haar measure ofUSp(2n,C). Let 



Tr(HP) 



Tr(A-P)+l 



W{M) = ±^p ifj is odd, and W(M) = ll ^ )+1 if j is even. Then 



\W < x ) 



.'■() 



2 dx 



0{j/n). 



Proof. The result is trivial if Aj > 2n + 1, so assume that 4j < In + 1. We 
apply Theorem 12, II to the exchangeable pair (W, W') with a = -^-\ — > an d 
will take the limit £ — >• in each term (keeping j,n fixed). By part 1 of 
Lemma 12. 21 and Lemma [4. 51 the first term is 0(VJ l n )- By part 2 of Lemma 
12.21 and Lemma 14.71 the second term is 0(j/n). By the Cauchy-Schwarz 
inequality and Lemma 14,6} 

E\W' - W\ 3 < y/E(W - W) 2 E(W - Wy = 0(£ 3/2 ). 

Thus the third term in Theorem 12.11 tends to as t — > 0, and the result 
follows. □ 

5. The unitary group 

In this final section, we treat the unitary group U(n, C). We let p\ be 
as in IScction "31 and ISection 41 and define the real valued random variable 
W = Pj -Jr- . The main result of this section is a central limit theorem 

/2? 



for W, with error term 0(j/n). To begin, we recall the following moment 
computation from [5J. 

Lemma 5.1. Let M be Haar distributed on U(n,C). Let (oi,a2,--- , ajt) 
and (&i,--- , bj~) be vectors of non-negative integers. Let Z%, ■ ■ ■ ,Zk be inde- 
pendent standard normal random variables. Then for all n > J2i=i( a i + h), 



E 



Y[Tr(M j ) a i -Tr(Afi) 
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Rains [20| (see also [14J] ) determined how the Laplacian acts on power sum 
symmetric functions. We require his formulas only in the following cases. 



Lemma 5.2. (1) 



(2) 



(3) 



&U(n)Pj = -njPj -J ^ PLj-l- 
l<Kj 



A U(n)Pj,j = -ZnJPjj ~ 2J 2 P2j ~ 2jPj Yl PU-l- 



Au{n) {PjPj) = 2 fn - ZnjpjPj - jpj Y pij_i - jpj Y Pl,j-l- 

l<Kj l<l<j 

Lemma 15.31 computes the conditional expectation E[W'|M]. 

Lemma 5.3. 

E[W'|M] = (1 - njt) W + R(M), 
with 



R{M) = t 



§ Y pi,j-i-]l J 2 Y pij-t 



+ 0(t 2 



Proof. Applying Lemma 12.31 and part 1 of Lemma 15.21 gives that 

E[W'\M] 
= e tA (W) 



W + t 



-njW 



i<Kj 



Y PiJ-i-Vo ^ Pl *- 1 



l<Kj 



+ 0(t 2 ), 



as desired. 



D 



Lemma 15.41 computes E[(W' — W) 2 |M]. As in the other cases, there are 
nice cancelations. 

Lemma 5.4. 

E[(W - W) 2 \M) = tj (2n - P2j ~m) + °( t2 )- 
Proof. Clearly 

E[(W' - Wf\M] = E[(W') 2 \M] - 2WE[W'\M] + W 2 . 
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By Lemmas 12.31 and 15.21 

E[(W') 2 \M] 

= W 2 + — A\p jtj + 2 Pj p~ + pj]} + 0(t 2 ) 



2j 



W 2 + t 



-npjj - jp2j - Pj Y PU-i ~ n PhJ ~ JPij 

l<Kj 



+t 



-pj Y py-' + 2 i n ~ 2n Pm - pj Y Pl >j- 1 ~ p i Y PtJ-t 

l<l<j l<Kj 1<1<3 

+0(t 2 ). 
By Lemma E3 -2WE\W'\M\ is equal to 



-2W 2 + t 



npjj + 2np jPj + n Pj j + Pj Y PU-l 

1<1<3 



+t 



Pj Y P l >i~ l + Pi Y Pl >i~ l + p i Y P^i- 1 

l<l<j 1<K? l<Kj 



+ o(t 2 



Thus 



E[(W') 2 \M] - 2WE[W'\M] + W 2 = tj [2n - p 2j -P2j] + 0(t 2 ) 
and the lemma is proved. 
Lemma 5.5. Suppose that Aj < n. Then 

Var(E[(W' - W) 2 \M]) = Aft 2 + 0{t 3 ). 
Proof. By Lemmas 15.41 and 15.14 

Var(E[(W' -W) 2 \M}) = j 2 t 2 Var(p 2j +p^) + 0(t 3 ) 

= j 2 t 2 E[(p 2j +p^) 2 ] + 0(t 3 ) 
= Aj 3 t 2 + 0(t 3 ). 



□ 



D 



Lemma 5.6. Suppose that Aj <n. 

(1) E(W - W) 2 = t2jn + 0(t 2 ). 

(2) E(W'-W) i = 0(t 2 ). 

Proof. Lemma El implies that E(W -W) 2 = E [tj (2n - p 2j - pTj)] + 0{t 2 ). 
From Lemma [5TTT E(p 2 j) = E(p 2 ~j) = 0, which proves part 1 (using only that 
2j<n). 

For part 2, first note that since 

E[{W'-Wf] = E{W 4 )-AE{W 3 W')+QE[W 2 {W') 2 ]-AE[W{W') 3 ]+E[(W')% 

exchangeability of (W, W) gives that 

E(W' -W) A = 2E(W 4 )-8E(W 3 W) + QE[W 2 (W') 2 ] 

= 2E(W 4 ) - 8E[W 3 E[W'\M\] + 6E[W 2 E[(W') 2 | M ]]- 
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Using Lemmas 15.21 and 15.31 this simplifies to 

2E(W 4 ) - 8E[W 4 } + 6E[W 4 } 

8njW* + 8W 3 J J - £ PU-i + SW 3 J J - J^ KFi 



+ffi 



l<Kj 



l<Kj 



+tE 



+tE 



+£E 



-6nW 2 pjj - 6jW 2 p 2j - 6W 2 pj ^ py_j - 6nW 2 pJJ 

1<1<3 



SjW 2 Wj - SW 2 p] J2 KFl + 12njW 2 



-12nW 2 Pj p- - 6W 2 P] Y, KJ~i-QW 2 pJ Y PlJ-l 
l<l<j l<Kj 



+ 0(t 2 ). 



By Lemma 15.11 after dropping out terms with expectation, there remains 

tE^W^jn - 6W 2 n PjJ - 6W 2 np-J- + 12W 2 jn - l2W 2 npjPj] + 0(t 2 ) 
= t[24jn - 6jn - 6jn + 12jn - 24jn] + 0(t 2 ) 
= 0(t 2 ), 
as needed. □ 



Lemma 5.7. Let R = t 



2 /^l<l<j PlJ-l Y 2 2-il<l<j Plj-l 



+ 0(t 2 ), 



and suppose that 4j < n. Then E[R 2 ] < 2-i — \- 0(t 3 ). 



Proof. Applying Lemma f5.1l and keeping only terms with non-0 contribution, 
one has that 

E[i,> 2 ] = jt 2 E[ J2 PiJ-iKFl] + 0(t 3 ). 



l<Kj 



If j is odd, then by Lemma loTTl 



E[R 2 ]=jt 2 ^[/(j-0] + O(t 3 ) 
while if j is even, one obtains that 



s,_U 4 -f) 



t 2 + o(t 3 



E[R 



, 2l _ (2j 4 + 3j 3 - 2f 
12 



-t 2 + 0(t 3 ). 



The result follows. 



□ 



W{M) = ^=[Tr{Mi)+Tr{Mi)}. Then 



Theorem 5.8. Let M be chosen from, the Haar measure ofU(n,C), and let 
\W < xq) ^= / e'^dx = 0(j/n). 
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Proof. The result is trivial if 4j > n, so assume that 4j < n. We apply 
Theorem 12.11 to the exchangeable pair (W, W') with a = tnj, and will take 
the limit t — > in each term. By part 1 of Lemma 12.21 and Lemma 15.51 
the first term is at most ^ . By Lemma 15.71 and part 2 of Lemma 12.21 
the second term in Theorem 12.11 is at most -^-. By the Cauchy-Schwarz 
inequality and Lemma EH 

e|w' - w\ 3 < ^/e(w - w) 2 e(w - wy = 0(t 3/2 ). 

Thus the third term in Theorem 12.11 tends to as t — > 0, and the result 
follows since 

12VJ , 19j < 22j 



n 2n n 
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